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Abstract 

New Ekpyrotic Cosmology is an alternative scenario of early universe cosmology 
in which the universe existed before the big bang. The simplest model relies on two 
scalar fields, whose entropy perturbation leads to a scale-invariant spectrum of density 
fluctuations. The ekpyrotic solution has a tachyonic instability along the entropy field 
direction which, a priori, appears to require fine-tuning of the initial conditions. In this 
paper, we show that these can be achieved naturally by adding a small positive mass 
term for the tachyonic field and coupling to light fermions. Then, for a wide range of 
initial conditions, the tachyonic field gets stabilized with the appropriate values well 
before the onset of the ekpyrotic phase. Furthermore, we show that ekpyrotic theory 
is successful in solving the flatness, horizon and homogeneity problems of standard big 
bang cosmology. Motivated by the analysis of the tachyonic instability, we propose 
a simplification of the model in terms of new field variables. Instead of requiring 
two exponential potentials, one for each scalar field, it suffices to consider a single 
nearly exponential potential for one of the fields and a tachyonic mass term along the 
orthogonal direction in field space. All other terms in the potential are essentially 
arbitrary. This greatly widens the class of ekpyrotic potentials and allows substantial 
freedom in determining the spectral index and possible non-Gaussianity. We present a 
generalized expression for the spectral index which is easily consistent with the observed 
red tilt. We also argue that for a wide range of potentials non-Gaussianity can be 
substantial, within the reach of current observations. 



1 Introduction 



The ekpyrotic scenario PQ [21 El S] is an alternative theory to the standard inflationary big 
bang paradigm. Instead of invoking a flash of exponential expansion shortly after the big 
bang, ekpyrotic theory proposes that the universe started out in a cold state followed by 
a long period of slow contraction. Despite this stark difference in dynamics, remarkably 
the two models make identical predictions for the origin of structure formation: a nearly 
scale-invariant, adiabatic and Gaussian spectrum of density perturbations. One testable 
distinctive prediction is that inflation also generates scale-invariant primordial gravitational 
waves whose amplitude in the simplest models is within reach of near-future microwave 
background polarization experiments whereas the ekpyrotic scenario does not [H E]. 
The idea of a pre-big bang contracting phase arising from nearly vacuous initial conditions 
originated in the pre-big bang scenario of Gasperini and Veneziano [7j. 

Until recently, ekpyrotic theory suffered from two important drawbacks. 

First, it was unknown how to describe a non-singular bounce within a consistent effec- 
tive theory without introducing ghosts or other catastrophic instabilities [8]. This lead the 
advocates of the Big Bang/Big Crunch ekpyrotic [21 [3] and cyclic models P [10] to propose 
that the universe undergoes a big crunch singularity to be resolved by stringy physics. This 
seems plausible in light of the mildness of the singularity. In the context of colliding branes 
in heterotic M-theory [TTl [T2] , for instance, the singularity corresponds to the fifth dimen- 
sion shrinking to zero size, with the large three spatial dimensions remaining finite. Despite 
intense activity in recent years [T3j, however, a definitive example of a cosmological bounce 
in string theory is still missing. 

Second, the fate of perturbations through the bounce, and therefore the prediction of scale 
invariance, is ambiguous. This stems from the fact that the scale-invariant growing mode of 
the scalar field fluctuations precisely projects out of C [31 HI] — the curvature perturbation 
on uniform-density hypersurfaces [15]. The latter is a useful variable to track since it is 
conserved on super-horizon scales. For this reason, many have argued that the outcome 
of the bounce would amount to matching ( continuously, resulting in an unacceptably blue 
spectral tilt [16J. While this is universally true for any non-singular bounce within 4d Einstein 
gravity, it is nevertheless conceivable that stringy or higher- dimensional effects relevant near 
the bounce may lead to mode-mixing and impart ( with a scale- invariant piece [17j. 

In a recent paper [TH], we proposed a fully consistent and complete scenario which ad- 
dresses both issues. While sharing many important ingredients of the old ekpyrotic sce- 
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nario [H Ellin], this model of New Ekpyrotic Cosmology generates a non-singular bounce, all 
describable within a 4d effective field theory, and unambiguously leads to a scale-invariant 
curvature perturbation. 

The physics of the bounce exploits the mechanism of ghost condensation [20] , which relies 
on higher-derivative corrections to the kinetic term for a scalar field, akin to k-inflation [21] or 
k-essence [22]. Ghost condensation allows for violations of the null-energy condition (NEC) 
without introducing ghosts or other pathologies [23]. In [IH] we showed explicitly how the 
ghost condensate can be merged consistently with the preceding ekpyrotic phase to generate 
a smooth transition from contraction to expansion. Whether ghost condensation can be 
realized in a UV complete theory, such as string theory, remains an open issue [2l] , although 
see [23] for a recent attempt. Even if a smooth bounce is realized by very different physics, 
many of our key results will still apply. For our purposes the ghost condensate provides an 
explicit and tractable realization of such a bounce. 

In our model ( acquires a scale-invariant spectrum well-before the bounce through the 
conversion of entropy (or isocurvature) perturbations. By having two scalar fields, each with 
their own ekpyrotic potentials, the entropy perturbation — corresponding to the difference 
in the scalar field fluctuations — is scale- invariant [IHl [2S1 123 I2H]- See also [23, EHl EI]- 
In [18], we showed how the entropy mode gets converted into the adiabatic mode by using 
features in the potential which are independently necessary for ending the ekpyrotic phase 
and bridging into the bounce. Then, since the physics of the bounce is all within 4d effective 
theory, ( goes through the bounce unscathed and emerges in the hot big bang phase with a 
nearly scale-invariant spectrum. 

As pointed out in [IHl [2S] and studied extensively in [321 E3], the model has a tachyonic 
instability along the entropy direction. Since this instability is at the origin of the growth 
of entropy modes, the tachyon is unavoidable [34]. As it stands. New Ekpyrotic Cosmology 
would ostensibly require fine-tuned initial conditions in order for the field trajectory to start 
out rolling along this tachyonic ridge. In this paper, we show how such initial conditions can 
be achieved naturally with a pre-ekpyrotic stabilization phase. 

Our original scenario [18J was cast in terms of two scalar fields, (pi and 02, each with their 
own steep, negative exponential potential. Through a rotation in field space, following [32j 
we can instead study the dynamics in terms of new field variables and respectively the 
field directions along and orthogonal to the field trajectory. This is reviewed in Sec. [2j In 
this language, the solution describes rolling down an ekpyrotic potential along 0, while x 
remains fixed at a tachyonic point, thereby making the instability manifest. 
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The (0, x) perspective leads us to propose a significant simplification and generalization 
of the scenario, described in Sec. |3| By treating and x as fundamental fields, the origin of 
the potential in terms of the two steep and quasi-exponential potentials becomes unnecessary. 
All we need is a single steep and nearly exponential potential along 0, as well as a tachyonic 
mass for x, leaving tremendous freedom in specifying the global shape of the potential. The 
only constraint which is crucial in generating a scale-invariant spectrum is that the curvature 
of the potential along the x 4> directions must be nearly the same: V^^^ ~ ^ax- 




Figure 1: General shape of the potential during the ekpyrotic phase. The arrow indicates 
the desired solution, corresponding to rolling down a steep, negative and quasi-exponential 
potential along 0, while remaining perched on top of a tachyonic ridge along x- 



This novel framework greatly expands the class of allowed ekpyrotic potential — see 
Fig. [T] for a generic example. In Sec. 3.2| we derive the spectral tilt for the most general 
New Ekpyrotic potential and find that it depends on three parameters: the usual fast- 
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roll parameters e and f] characterizing respectively the steepness and deviation from pure 
exponential form of the potential, as well as a new parameter 6 describing the difference in 
curvature between the x and (p directions. The class of potentials studied in [18] corresponds 
to 5 = 0. Here we see that turning on 6 allows for even greater freedom in the spectral 
tilt. In particular, pure exponential potentials [i] = 0) were found to give a blue spectrum 
in [m Eg ETj, in disagreement with the recent WMAP data [35j. This led [T8| 126] to 
consider deviations from the pure exponential form (77 7^ 0), which does allow for a small 
red tilt. Here, however, we see that even for pure exponentials the tilt can be red by having 
non- vanishing S. 



In Sec. 3.3, we discuss the generic form of higher-order terms in x- These self-interaction 
terms play a crucial role in determining the level of non-Gaussianity of the perturbation 
spectrum. We find that the non-Gaussian signal generated during the ekpyrotic phase is 
generically substantial, within the reach of current and near-future experiments. 

The tachyonic instability in the x direction is part of a more general discussion of initial 
conditions in New Ekpyrotic Cosmology, in particular the assumed degree of initial flatness 
and homogeneity. In the context of the original ekpyrotic scenario [T] it was initially believed 
that the model did not address the homogeneity and flatness problems of standard big bang 
cosmology [36] • See [37] for a recent critique of the scenario along these lines. Such criticism 
motivated the advent of the cyclic model [9l [10], where homogeneity and flatness in a given 
cycle is achieved through a phase of dark energy domination in the previous cycle. It was 
recently realized, however, that the ekpyrotic model does in fact solve these problems [38j . 
in a way akin to inflationary cosmology. 

In Sec. |4| we argue that New Ekpyrotic Cosmology addresses the flatness, homogene- 
ity and isotropy problems of standard big bang cosmology without invoking any phase of 
cosmic acceleration. This relies crucially on the ekpyrotic attractor mechanism — during 
the ekpyrotic phase, the energy density in the ekpyrotic scalar field blueshifts much faster 
than any other component, in particular spatial curvature and anisotropic stress. In other 
words, contrary to the naive expectation that a contracting phase makes the universe more 
inhomogeneous, here the universe becomes increasingly smooth. Therefore, making the nat- 
ural assumptions that the proto-ekpyrotic patch is fairly homogeneous, isotropic, and has 
curvature comparable to the Hubble scale, then, provided the ekpyrotic phase lasts for suf- 
ficient number of e-folds, the universe emerges in the hot big bang phase with the high level 
of symmetry we observe today. As such. New Ekpyrotic Cosmology resolves the standard 
problems of the big bang model. 
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This is precisely analogous to inflationary cosmology. There one assumes that over some 
macroscopic patch the universe is sufficiently homogeneous and ffat to allow inffation to take 
over. Subsequently, cosmic acceleration ffattens and homogeneizes the universe to the high 
level we observe it today. Of course, in absolute terms the assumed initial radius of curvature 
is much larger in ekpyrosis than in inffation, but by no means is this a drawback of one model 
versus the other. Instead, this corresponds to drastically different assumptions about the 
initial state of the universe: inffation assumes a hot initial state, with correspondingly high 
Hubble scale; ekpyrosis proposes a cold beginning, with correspondingly small expansion 
rate. In the absence of a concrete and rigorous theory of initial conditions, either starting 
point is equally reasonable. With respect to ffatness and homogeneity, the proper yardstick 
is, therefore, the assumed level of ffatness and homogeneity compared to the respective initial 
Hubble radius. And, from this point of view, both models are equally successful. 




Figure 2: A pre-ekpyrotic mechanism sets up the desired initial conditions for the ekpyrotic 
phase. With the addition of a positive mass-squared term, the ffeld is initially stable in the x 
direction. By introducing couplings to light fermions, the ffeld settles down to the minimum 
for a wide range of initial conditions by the onset of the ekpyrotic phase. 

In this paper, we address in a simple and natural way the ffne-tuning of initial conditions 
associated with the aforementioned tachyonic instability of the two-ffeld model. In Sec. |5]we 
quantify the nature of the instability and discover that exponentially ffne-tuned conditions 
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must be satisfied at the onset of the ekpyrotic phase in order for the roll along the tachyonic 
ridge to last sufficiently long to produce enough e-foldings of scale-invariant perturbations. 
We then solve this initial fine-tuning by introducing a positive mass squared term for x, which 
is relevant and stabilizes x &t early, pre-ekpyrotic times, but becomes negligible during the 
ekpyrotic phase. In other words, this term serves to set up the desired initial conditions but 
does not jeopardize the subsequent generation of perturbations. Thus, x starts oscillating 
around its stable minimum. By introducing couplings to light fermions, the energy in these 
oscillations is quickly converted into thermal radiation. We find that for a wide range of initial 
conditions, spanning many orders of magnitude in initial x energy density, x is exponentially 
close to the top of the tachyonic ridge by the onset of the ekpyrotic phase. The modified 
potential is sketched in Fig. [2] 




Figure 3: The end of the ekpyrotic phase is triggered by a term in the potential which is 
small at early times but eventually pushes the field away from the tachyonic ridge. The field 
rolls down to a minimum in the x direction and starts oscillating around it. Meanwhile the 
field keeps rolling along the direction. 

In Sec. |6} we turn to a discussion of the exit from the ekpyrotic phase, the subsequent 
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NEC- violating phase and the resulting non-singular bounce. In this part of the story as well, 
the {(f), x) picture leads to tremendous simplications compared to our original model [18] . 
To bring the ekpyrotic phase to an end, we include corrections terms which eventually drive 
X away from the tachyonic ridge and towards a stable minimum, as shown in Fig. |3| This 
is necessary to generate a turn in the field trajectory, thereby imprinting the scale-invariant 
entropy perturbation spectrum onto (. 

Thus X ^^oUs to the minimum and starts oscillating around it. Its mass around the 
minimum is generically large compared to Hubble, and therefore the energy density in x 
blueshifts as a dust component. Note that the oscillations are not damped by the light 
fermions of the pre-ekpyrotic phase since, as x rolls away from the tachyonic ridge, these 
fermions become heavy and can be henceforth ignored. Nevertheless the energy density 
in X quickly becomes subdominant as (f) keeps on rolling, due to the ekpyrotic attractor 
mechanism. The dynamics therefore effectively reduce to single-field ekpyrotic motion along 

The rest of the story consists of entering a ghost condensate phase which leads to 
a violation of the NEC and results in a non-singular cosmological bounce. As described 
in [IH], this requires the post-ekpyrotic potential to become positive and fiat. Thus, shortly 
after x rolls off, we envision reaching a minimum of the potential followed by a steep 
rise to positive values where the potential becomes fiat. (Note that although the potential 
is fiat and positive, there is no sustained inflationary phase here. From the onset of the 
ghost condensate phase until reheating, the scale factor changes only by a factor of order 
unity.) Climbing up to this plateau greatly reduces the kinetic energy in 0, which allows 
to enter a ghost condensate phase. While most of this bounce story parallels that of our 
original scenario [18], a key difference here is that we only need one ghost condensate field, 
as opposed to two. This is another considerable simplification brought about by the 0, x 
framework. 

The successful merger of ekpyrosis and ghost condensation leads to a number of con- 



sistency conditions which we derive explicitly in Sec. 6^ and summarize in Sec. 6^ For 
example, the energy density in the x oscillations must be small enough to allow the ghost 
condensate to dominate the dynamics and lead to a bounce. None of these conditions are 
very constraining, however, and we illustrate how they can be satisfied for a wide range of 
parameters by studying a specific potential in Sec. [7} 

In the Conclusion in Sec. |8| we step back and assess the current status of New Ekpyrotic 
Cosmology compared with the inflationary model, both at the theoretical and phenomeno- 
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logical levels. 

The sequence of events in New Ekpyrotic Cosmology is outlined in Fig. |4j Starting 
from a cold initial state at some initial time ti, the field x begins oscillating around the 
stable minimum and decays into radiation — see Fig. |2] Since radiation blueshifts more 
slowly than the energy density in 0, eventually the latter takes over and dominates the 
energy. Meanwhile, the x direction becomes tachyonic, allowing for the growth of entropy 
perturbations. This marks the onset of the ekpyrotic phase, denoted by tek-beg- During this 
phase, the universe contracts very slowly, so that a{t) is nearly constant, whereas the Hubble 
radius shrinks by an exponential amount. It is also during this phase that fluctuations 
in X are amplifled and acquire a scale-invariant spectrum. The ekpyrotic phase comes to 
an end at tek-end, when x is pushed away from the tachyonic ridge and starts oscillating 
around a new minimum — see Fig. |3] Shortly thereafter {t = tc), 4> reaches a minimum in 
the potential, climbs up on the flat positive plateau, and enters the ghost condensate phase. 
This violates the NEC and leads to a non-singular bounce. Not long after the universe starts 
expanding, the energy density in is converted into matter and radiation, which reheats 
the universe and marks the beginning of the hot big bang phase {t = treheat)- The rest of 
the story until today (t = to) is standard: the universe successively undergoes epochs of 
radiation, matter and dark energy domination. 

stabilization ekpyrotic x moves NEC-violating phase hot big bang phase 

phase phase to new min and bounce 



^ek-beg 



Scale factor 



Approx. const. Approx. const. Changes by factor Grows by 



a.( t) of order unity 



Hubble radius ^ Shrinks by e^^k Approx. const. Grows infinitely big Grows by e^^rad 

// ^t) but comes back to original 



Figure 4: Sequence of events in New Ekpyrotic Cosmology. 



2 Two-Field Ekpyrosis 

The ekpyrotic mechanism for generating a pre-big bang, scale-invariant spectrum for (, as 
described in [TH], relies on two scalar flelds, 0i and 02, each rolling down a steep, negative 
and nearly exponential potential. For concreteness, let us assume for now that the potentials 
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are pure exponentials, in general with different powers pi <C 1 and p2 "C 1: 

n*>.)^-f.exp(-yi^)-V,exp(-yi^) . (2.1, 

In this paper, Mp\ is chosen to be the "reduced" Planck mass where Mpi = 2.4 x 10^^ GeV. 
Potentials that deviate from the pure exponential form were discussed in [iHl [26]. In this 
section, the two scalars can be taken to have canonical kinetic terms; all higher-derivative 
self-interactions, necessary to generate a ghost condensate and produce a smooth bounce, 
are assumed to be negligible during this ekpyrotic phase. Under these assumptions, 0i and 
02 satisfy the standard scalar equations of motion 

01 + 3if0i = -V^, ; 

02 + 3if02 = -V^, . (2.2) 
These equations, combined with the Friedmann equation, allow for a scaling solution: 

a{t) ~ {-ty^+P^ ; H = Pl±Pl ■ 



01 (t) = v^Mpilog 




Vi 



M2jpi(l-3(pi+p2)) 



02(t) ^ v/2p.Mp.log I -J m^Mi\.^P2)) ' ' ' ^'-'^ 



where t is chosen to be negative and increasing as t — > 0. For pi,p2 <^ 1, (2.3) describes a 



slowly-contracting universe driven by a highly stiff fluid with equation of state 

w = —^ r-l>l. (2.4) 

3(Pl+P2) 

This cosmology is the essential feature of ekpyrotic dynamics and is the origin of the gener- 
ation of a scale-invariant spectrum in this context. 

Ekpyrotic dynamics was first studied for a single scalar field, in which case the scaling 
solution has the property of being an attractor. Indeed, since w ^ 1 for the background, 
the corresponding energy density blueshifts much faster than all other relevant components 
— curvature, matter, radiation, coherent kinetic energy and anisotropic stress. This is the 
precise analogue in a contracting universe of why accelerated expansion is an attractor in an 
expanding universe. In the latter context, the nearly constant potential energy of a slowly- 
rolling scalar field comes to dominate because it redshifts slower than any of the above 
components. 
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2.1 Tachyonic Instability 



In the two-field case, however, the scahng solution (2.3) is not an attractor because of 
a tachyonic instability in the direction orthogonal to the field trajectory — the so-called 
entropy direction. This was first pointed out in [THl [2S] and studied extensively in [321 
33]. This instability is essential to the generation of a scale- invariant spectrum of entropy 
perturbations and, hence, cannot be removed [iHl [31]. However, as we show later in this 
paper, it is easy to include further terms in the potential which are negligible during the 
ekpyrotic phase, but which at early times bring the field trajectory exponentially close to 
the desired one. 

The tachyonic instability is most easily seen in terms of two new field variables, (j) and 
X, given by 



X=- ; , , (2.5) 



The scaling solution above implies that x remains constant on the background trajectory. 
Denoting this constant by Xti we find that 



At Xt the potential is tachyonic in the x direction. To see this, note that the potential (2.1) 
can be rewritten in terms of and x as 



V^(0, X) = 'Voe-r^ "'"^^ ( ^e'/f VK--)/^-^" + ^^Vf VK--)/^" ) , (2.7) 



where p = pi + p2 and Vq = (1 + ^2/^1)^1 exp {^^/jhfihP Xt/M-pi^ . Taylor-expanding (2.7) 
as a power series in x — x^, we find to quadratic order that 



V{<P.x) = -Voe-^-^'*'"'^'[l + ^{x-Xt? + ...j , (2.8) 

thereby revealing a tachyonic mass squared in the x direction. This tachyonic ridge is 
sketched in Fig. [T] Note that (2.8) satisfies 



Vxx = Vh (2-9) 

at X = Xti an essential feature in generating a scale-invariant spectrum for 5xi as we review 
below. 



In terms of and x, the scalar equations of motion (2.2) take the form 

+ 3i/0 = -V4, ■ 

X + 3Hx = -V^, (2.10) 
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while the background solution (2.3) becomes 

ait)- 



-ty 



H 



p 
t 



X = Xt 



2pMpilog 



Vn 



(2.11) 



This describes the system rolling along the direction on the top of the tachyonic ridge in 
X with the effective scalar potential for given by 
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Kfr(0) = -V'oexp 



PI 



(2.12) 



Note, using (2.8), (2.11) and (2.12), that 



V, 



(2.13) 



_ 2 

' ,<P<P\X=Xt ~ ^ff,</><^ ^ 

for p <^ 1. As rolls down this steep, negative exponential potential, the tachyonic mass 
for X g6ts larger in magnitude. A useful relation between the latter and other background 
quantities is 



2 =_T/ I ^ ^ 

"^tachyon = ~Kxxlx=Xt ~ 



(2.14) 



where we have used (2.9), (2.11) and (2.13). Since p <C 1, the rate of instability is always 



much greater than the Hubble parameter. 



2.2 Instability and Scale-Invariant Entropy Perturbation 

Since the field trajectory is orthogonal to X; the fluctuations in the latter by definition 
coincide during the ekpyrotic phase with the entropy perturbations. Moreover the tachyonic 
mass uncovered above is essential in amplifying these perturbations, as we now review. This 
implies that the instability cannot be cured without spoiling the scale invariance of the 
entropy perturbation spectrum. 

The evolution of the Fourier modes of 6x is governed by 



Sxk + SH5xk + {K + V,^^]Sx 



0. 



(2.15) 



In Sec. 3.2 we will solve (2.15) and derive a general expression for the spectral tilt, including 



Hubble damping and corrections from the potential. For the moment, we note from (2.14) 
that 

P 

V2 



H 



(2.16) 



XXI 
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at X = Xt- Thus, since p -C 1, we can neglect the Hubble damping term in (2.15) and treat 
space as fiat — we can set a = 1 without loss of generality. In this approximation, x is a 
free scalar with time- dependent mass: 







(2.17) 



where we have used (2.9) and (2.13) to replace V^^ by — 2/t . 



Well within the horizon, a given fc-mode starts out in the usual Bunch-Davies vacuum: 



5xk = e jylk. The solution to (2.17) with these initial conditions is 



-ikt 



^Xk 



1 



ki 



(2.18) 



On super-Hubble scales, k{—t) ^ 1, up to an irrelevant phase factor this reduces to 



V2i-t) 



(2.19) 



corresponding to a Harrison-Zeldovich spectrum. It follows that, as stated above, the scale 
invariance of the entropy spectrum in two-field ekpyrosis relies on the existence of a tachyonic 
instability along one field direction, denoted by x above. 



3 Generalized Ekpyrotic Potentials 

From the above analysis we learn two important lessons about the ekpyrotic phase: first, the 
physics associated with two exponential potentials for (pi and 02 is most transparent when 
the theory is rewritten in terms of the rotated fields and second, that a nearly scale- 
invariant spectrum arises in the fluctuations of x as a result of the relations V^^ ~ 
and V^^ = V^^ at x = X* potential (2.8). From the first fact we conclude that it is far 



simpler to treat and x as the fundamental fields in the theory, which we do henceforth. It 
follows from the second fact that we are free to modify the potential in the ekpyrotic phase 
as long as the relations 

V^^^-2/t\ K00^Kxx (3-1) 



at X = Xi continue to hold. This new perspective greatly simplifies New Ekpyrotic Cosmology 
and broadens the class of allowed potentials. 
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3.1 Generalized Potentials 



The simplest such modification is the following. Motivated by (2.8), consider any potential 
of the form 



V^(0,x) = -^oe V? 



1 + 



1 



ix-Xt? + Fi<P,x-Xt) 



(3.2) 



where Xt is any constant. Meanwhile, the function F satisfies 



,XX\X=Xt 



(3.3) 



but is otherwise arbitrary. Since the exponential in and the quadratic term in the Taylor 



expansion in x — Xt are identical to those in (2.8), this modified potential clearly continues 



to satisfy conditions (3.1). Note that in the case of two exact exponentials discussed above. 



keeping higher order terms in the expansion (2.8) gives 



Fi(j),X- Xt) 



(3.4) 



which manifestly satisfies (3.3). Hence, (2.8) is indeed of the form (3.3). Since the generalized 



functions F encode the self-interactions for x, they play an important role in the analysis of 
non-Gaussianity in New Ekpyrotic Cosmology, which we discuss below. 

The (0, x) perspective suggests an even broader class of allowed ekpyrotic potentials. 
Note that although the exponential in and the quadratic term in the Taylor expansion 
in X — Xt imply that the potential (3.2) satisfies conditions (3.1), they are hardly the only 



functions to do so. For simplicity of notation, let us henceforth set 



= 



(3.5) 



without loss of generality. Then, consider any potential of the form 



V(0,X) = V(0) ( 1 + ^^/(0)x' + i^(0,x) 



where V(0) satisfies the standard fast-roll conditions of single-field ekpyrosis 
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7] = 1 



PI \y 



V 



,<t> 
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< 1 



|?7| < 1 



(3.6) 



(3.7) 



which require the potential V(0) to be steep and nearly exponential, respectively. Meanwhile, 
the function /(0) has the property 

^ 1 (3.8) 



over the relevant range of modes, and F is constrained by condition (3.3). It is straightfor 



ward to show that potential (3.6) satisfies the conditions (3.1) for a nearly scale-invariant 
spectrum and is the most general potential to do so. Clearly the potential in (3.2) — with 
Xt set to zero — is of this form, with V(0) = -Voe~V^ ^/^-^p^ and = 1. 

One might object that, at first sight, the (0, x) perspective requires an unnatural fine- 
tuning — f{(f)) ^1 — while this is seemingly automatically satisfied in the original (f)i,4>2 
language. This would seem to suggest that the latter is a more convenient set of fundamental 
fields to study. On the contrary, we argue that the original picture in terms of (pi and 02 



requires much more fine-tuning. The generalized potential V{(j),x) in (3.6) requires one 
function to be nearly exponential and one function in the Taylor expansion — /(0) — to 
be nearly unity. However, the shape of the potential away from x = 0, given by the general 
function F, can be quite arbitrary. In contrast, in the 0i and 02 framework a second function 
must be of nearly exponential form, corresponding to tuning infinitely many coefficients in 
a Taylor expansion. 

The generalized function V(0) and the modified quadratic term in the Taylor expansion 
in X have a profound impact on the form of the spectral index of the theory. This will now 
be discussed in detail. 



3.2 Spectral Index for General Potentials 



In this section we derive the spectral tilt for the general potential (3.6 ). A necessary condition 
for scale invariance, as we will see, is that the function /(0) be approximately unity over the 
relevant range of modes, which ensures that V^^ ~ V^j,,/, near x = Xt = 0. Therefore, let us 
parametrize / as 

/(0) = l + 35; |5|«1, (3.9) 

where the factor of 3 is introduced for convenience. Note that, generically, S can be a function 
of 0. For comparison, the analysis of Sec. 4 in |THj focused on the special case of identical 
potentials for 0i and 02, which translates in the 0, x variables to /(0) = 1. Turning on a 
non-zero 6 is not only more general but, as we will see, allows for greater freedom in the 
spectral tilt. 
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The perturbation equation for the entropy field 5x is given in (2.15) by 

6xk + SH6xk + + Kxx) SXk = 0. (3.10) 



When (5 = 0, we find from (3.6) that = at x = 0, in which case the calculation is 
identical to that presented in [TSl [20] • When 5 7^ 0, however, V^^ = V(j,(j,{l + 35) at x = 
and thus ^ 

6xk + ^H6xk + + + 35)) Sxk = 0. (3.11) 

To calculate n^, we follow the approach developed in [T^. Define the equation of state 
parameter 

3, , H d\nH 
e = -(l + w) = - = , 3.12 



where = In a is a dimensionless time variable. Equation (3.11) can now be recast entirely 
in terms of e and its derivatives with respect to A^, as shown in [18] for the 5 = case. Since 
the calculation closely parallels that of [18j, for the sake of brievity we only outline the key 
steps and refer the reader to [H] for details. 

In the ekpyrotic phase, the equation of state parameter is large: e ^ 1. Moreover, for 
simplicity, we assume that it is slowly-varying in time. In this limit, the e and t] fast-roll 
parameters can be written as 

1 ^ de 

Terms involving higher-derivatives on e, such as {de/dNy and d'^e/dN'^, are found to be 
higher-order in e and r] and, therefore, can be consistently neglected in the fast-roll approx- 



imation. To illustrate how this substitution works, we note, for example, that V^f,^ in (3.11) 
can be rewritten as 

+ + + + f^-") 

where the ellipses stand for terms of order e^^Tf^et] and so on. Let us emphasize that e 
and 7] can therefore be treated as essentially constant in this approximation, as in slow-roll 
infiation. Their time-dependence generate order e^^yf^er] corrections, which are negligible. 
The next step consists of introducing a new time variable 

and rescaling 5xk as 

Vk = a5xk- (3.16) 
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After some algebra, the perturbation equation (3.11) can be cast in the form 



+ ^^vu - 2 (1 - 3(6 - - 5)) = , (3.17) 

where we have dropped higher-order terms in the fast-roll parameters. 

Thus far, e and r] are approximated as constant whereas 5 is up to this point an arbitrary 
function of 0, and thus an arbitrary function of time. In order for the spectrum to be 
nearly scale- invariant, however, clearly 5 cannot vary rapidly while the observable modes 
exit the horizon. A small time-dependence is allowed, of course, which would lead to a 
running spectral index. We leave this interesting possibility aside and henceforth focus on 
the simplest case where 5 is approximately constant over the relevant spectral range. 

In terms of ?/ = s^/Vl "~ 8//, the solution with Bunch-Davies vacuum is given as usual by 
a Hankel function: 

''^ = ^{\^^n\y). (3.18) 



with n = (3/2)a/1 - 8(e - - S)/3 ^ 3/2 - 2(e - r/ - 5). Note that in the limit of fiat space 



(a ^ 1) and exact scale-invar iance {n 3/2), this reduces to (2.18). 

On large scales, k <^ aH, the amplitude tends to Vk ~ k~"-, corresponding to a spectral 
index: 

-1 = 4(e-r7-5). (3.19) 

As a check, for 6 = this reduces to the result of [TS]. Thus, 6 acts effectively as a correction 
to the fast-roll parameters contributing to the spectral tilt. 

As discussed in [TSl [23, 122], pure exponential potentials — 77 = 5 = in our language — 
yield a blue spectrum. However, as shown in [TSl ES] for the case 5 = 0, it is easy to get a red 
tilt by considering potentials that deviate from the pure exponential form, corresponding to 
non-zero 1]. For example, V{(f)) ~ exp{— </)"'), with n > 2, gives a red tilt at large 0. What we 
have found here is that the more generic case with 5 7^ allows greater freedom in getting 
a red tilt. Even if e — r/ is positive, for instance, which by itself would yield a blue spectral 
index, it is still possible to get a red tilt if 6 is positive and dominant. 

3.3 Non-Gaussianity in General Potentials 



Having derived the spectral index for the generalized potential (3.6), we now turn to the 
question of non-Gaussianity of the fluctuations. Whereas the spectral index is determined 
by the dependent over-all factor V(0) and the form of the quadratic term in the Taylor 
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expansion in x, the level of non-Gaussianity is associated with the function F, which encodes 
all self-interaction terms for x- Here we provide an estimate of the non-Gaussianity level 
based on the generic size of these interaction terms. For simplicity, we shall ignore the effects 



of gravity, as we did in Sec. |2.2| to estimate the 5x 2-point function, and work to leading 
order in the fast-roll parameters e and rj. 

While it is clear from the above discussion that F is completely arbitrary, nevertheless we 
can give some naturalness arguments for the cubic and higher-interaction terms in x based 



on the form of the mass term. To leading order in e,?7, it follows from (3.7) that 

'^M _ -1 
V 



(3.20) 



Thus the quadratic term in (3.6) can be written as 

1. 



where A is an effective cut-off scale: 



A = Mpiv^. 



(3.21) 



(3.22) 



From this point of view, we expect higher-order interactions to be suppressed by the 
same scale, leading to a generic F of the form 

a{<P{t)) 



F{(t^.x) 



A3 



(3.23) 



Here the coefficient a is a priori an arbitrary function of 0. However, to simplify the analysis 
we shall take it to be a constant which, naturally, is expected to be of order unity. 

The three-point function for the fluctuations 5x is determined by the interaction Hamil- 
tonian [39], which we can read off from the cubic part of the potential: 



vm)) 



a 5x^ 
3 A^ 



a 



3At' 



Note that in the last step we have substituted (3.20) and used V, 



(3.24) 

a —2/t^. Since i^int is 

exactly of the form studied in [27], we can immediately apply their results to our calculation. 
We find that the 3-point function for the curvature perturbation C, at the end of the ekpyrotic 
phase {t = tek-end), inherited from 6x, is given by 



(Cfc^(tek-end)Cfc2(^ek-end)Cfc3(^ek-end)) ~ (27r)'^5 | kj 
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n,A;f 23/2A(-tek-cnd) 



(3.25) 



where ki denotes the magnitude of fcj, and all modes are assumed to be super-Hubble: 

3k ' 



^ok-end) <^ 1- Moreover, is the amplitude of the power spectrum, A^ = k^(l, which 



AJ/^ ^ 10"^ . (3.26) 



is fixed by WMAP 

Note that the shape of the 3-point function is of the local form |l0] . 

The level of non-Gaussianity is usually expressed in terms of the /nl parameter [ITl H2 
defined in terms of ( by 



C{x) = Ux)-p^Ugix), (3.27) 



where (g has a Gaussian spectrum. From (3.25) we immediately read off that 

_ 5^2 1/2^^^^ 5 1/2 a l^ck-cndl 



where we have substituted (3.22) and (2.11), taking into account that p 2e in our approx- 
imation. Note that the above sign ambiguity has to do with the details of how the entropy 
perturbation is converted onto C. [27], an issue we ignore for the present discussion. 

The last factor in (3.28) can be expressed in terms of the amplitude of density perturba- 
tions as [H] 

^1/2 ^ ^ ^heat ,3 29) 

where /3 is a model-dependent prefactor having to do with the details of the exit from the 
ekpyrotic phase. As we will see in Sec. |6| in the approximation that the ekpyrotic phase ends 
through a sharp turn in the field trajectory, then we have f3 = A6, where A^ is the change 
in angle in the trajectory in field space — see also (5.5) of [IBj. Thus (3 is at most of order 
unity in that case. 

Save for the (3 coefficient, this expression is identical to its counterpart in inflation, with 
e replacing the einf slow-roll parameter. Now, as indicated in Fig. |4] and as we will discuss 
later in detail, the magnitude of the Hubble parameter at the end of the ekpyrotic phase is 



roughly the same as at reheating: |ifek-end| ~ -f^reheat- Using this fact and substituting (3.29) 



into (3.28), we find 

/NL^T^^e-^ (3.30) 

At the 2cr level, the recent WMAP data constrains /nl to be within the range [351 133]: 
—36 < /nl < 100. Using the liberal end of this bound, we get 

e>2-10-3«_ ^33^) 
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Since e -C 1, the level of non-Gaussianity tends naturally to be large in New Ekpyrotic 
Cosmology. For example, we will see in Sec. 4.4 that p ~ 10"^ for Tj-eheat = 10^^ GeV, 



corresponding to e ^ 5 ■ 10~^. Furthermore, as discussed in Sec. |6| typically /3 ~ 
Hence, taking a ~ 0{1) yields /nl ~ 40, which is below, but not far from, the present upper 
bound. Lower reheating temperatures correspond to even smaller values of e and therefore to 
larger non-Gaussian signal. Of course, we are free to choose the parameter a to be small and, 
hence, ekpyrotic theory can always lie within the observational bound. The implications for 
current and future experiments will be discussed elsewhere. 

The significant non-Gaussian signal is a distinguishing feature of New Ekpyrosis compared 
to slow-roll inflation, as the latter is well-known to give unobservably small /nl- Significant 
non-Gaussianity in inflation is achieved in multi-field models in models with higher- 
derivative kinetic terms [2T1 |l5l HG], HT], as well as in scenarios where density perturbations 
are generated by another light field, such as the curvaton |3H] and modulon scenarios [49]. 
The large non-Gaussianity level of New Ekpyrotic Cosmology is also in sharp contrast with 
single-field ekpyrosis, where density perturbations are even more Gaussian that in slow-roll 
inflation |50]. 



4 Flatness and Homogeneity in Ekpyrotic Theory 

The tachyonic instability in x is required in order to generate a nearly scale-invariant fluc- 
tuation spectrum during the ekpyrotic phase. Its existence, however, necessitates a careful 
study of both the initial conditions of ekpyrotic cosmology, as well as the mechanism for 
ending the ekpyrotic phase and generating a Harrison-Zeldovich spectrum in the curvature 
perturbation. We will give a detailed analysis of both of these issues later in this paper. 
Here, however, we note that initial conditions are closely related to the required degree of 
flatness and homogeneity of the universe. In this section, we will show that both of these 
fundamental problems are naturally solved in New Ekpyrotic Cosmology. 

When the original ekpyrotic scenario was proposed, it was believed that the model failed 
to address the flatness and homogeneity problems of standard big bang cosmology. This was 
viewed as a drawback for ekpyrotic theory as compared to inflation [36l [37] . This concern 
motivated the cyclic extension of the scenario [HI [TO] , where flatness and homogeneity at the 
onset of each contracting phase is achieved by an epoch of dark energy domination at the 
conclusion of the previous cycle. However, much progress has been made since then and the 
time seems right to revisit these issues in the pure, non-cyclic ekpyrotic context. We will 
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argue that, contrary to initial beliefs, ekpyrotic theory is equally successful as inflation in 
addressing the standard problems of big bang cosmology [38j , without having to invoke any 
phase of cosmic acceleration. 



4.1 Contrast with Inflation 

Inflationary and ekpyrotic cosmology are based on starkly different assumptions about the 
initial state of the universe. 

In inflation, one envisions the universe as starting at the big bang in a hot, rapidly 
expanding and chaotic state. While most of the universe is wildly inhomogeneous, the 
assumption is that there exists a Hubble patch somewhere which is sufficiently smooth to 
allow inflation to occur. A short phase of cosmic acceleration then blows up this tiny region 
into a large, homogeneous, isotropic and flat universe. 

Ekpyrosis, on the other hand, proposes that the universe begins in a quiescent, nearly 
vacuous and cold state. In the context of the original ekpyrotic scenario, this starting 
point was motivated by nearly supersymmetric or BPS initial conditions. In this scenario, 
all physical scales at the onset of ekpyrosis — initial Hubble parameter, curvature of the 
universe, radiation temperature — are very small in absolute, particle physics terms. 

In the absence of a concrete theory of initial conditions, however, the question of whether 
the universe started out in a hot or cold state belongs to the realm of metaphysics. Thus, 
with our current understanding, a tiny initial Hubble parameter is no less natural than a 
GUT-scale Hubble. From a cosmology perspective, the relevant yardstick for assessing the 
degree of genericity of initial conditions should, therefore, not be based on the absolute 
scale of any observable, such as the initial spatial curvature, but rather on whether it is 
fine-tuned relative to other observables, such as the energy density in other components. We 
illustrate this point by discussing the flatness and homogeneity problems of standard big 
bang cosmology. 

4.2 Flatness and Homogeneity Problems in Inflation 

Viewed as an additional energy component in the Friedmann equation, the curvature of the 
universe contributes a fractional amount 



1 




a 
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Since both a and H evolve in time, evidently so does flk- In particular, from the onset of the 
hot big bang phase until the present time the scale factor grows by a factor of e^'^''. Now, 
assuming for simplicity that our universe has been radiation-dominated ever since reheating, 
we have if~l/t~l/a^. It follows that over the same period the Hubble parameter shrinks 
by e~^^"^«d In other words, aH changes by e"^'''''^. A useful expression for A^^ad is therefore 

7.7- 1 / C^reheat -^reheat \ 1 / ^reheat \ /, r,\ 

where the last step follows from the standard redshift relation a ~ 1/T. For instance, since 
the current temperature of the universe is Tq = 2.7K ~ 10^"^ eV, then an initial temperature 
of Treheat = 10^^ GcV yields the standard Nj-^d ~ 60, whereas a lower reheat temperature of 
10^ GeV corresponds to A^rad ~ 46. Incidentally, we immediately recognize e^'^^^'* as the ratio 
between the size of our observable Hubble patch in comoving units (Aq — Hg^/ao) to that 

at reheating (Areheat = -^reheat/preheat)- 

Therefore, during the radiation-dominated epoch flk grows by an exponential factor: 

q(0) 2 717-2 

■^'fe "reheat -'-'reheat „2N, 



^(reheat) ^2^2 



e'''^-'^ . (4.3) 



Since we know from observations that Qk is at most a few percent today, we are forced to 
conclude that Qk must have been exponentially small relative to all other components at the 
onset of the hot big bang phase. This is the flatness problem. 

Inflation addresses this fine-tuning through a short phase of cosmic acceleration. During 
this phase the scale factor grows by an exponential amount, whereas the Hubble parameter 
remains essentially constant. Thus we can define the total number of e-foldings of inflationary 
expansion from some initial time ti until reheating ^reheat as 

Tvr 1 / ^reheat -^reheat \ , / C^reheat \ /, i\ 

"^""■"( a,H. 

By definition e^'"*" measures the size of the proto-inflationary patch in comoving units (A, = 
H^^ /ai) to that at reheating. 

This immediately implies that any initial curvature is exponentially suppressed by the 
time of reheating: 

Q (reheat) /,2 7x2 

_ „-2Ari„f / . r\ 

^{i) - 2 TT2 ■ \^-'^) 

ii;;. "reheat-'-' reheat 

Then, as long as A'inf > Ajadj we find that ^^^^ < In other words, starting from the 

natural choice fl^^^ ~ ^(1); ^ long enough inflationary phase results in a sufficiently large 
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radius of curvature for the observable universe. This is the essence of how inflation solves the 
flatness problem. Instead of requiring an exponentially large initial radius of curvature, as in 
standard big bang cosmology, inflation tolerates an initial curvature component comparable 
to all other forms of energy density. 

The homogeneity problem follows automatically from the above considerations, pro- 
vided that space-time is essentially homogeneous within the proto-inflationary patch. If 
Nini ^ A^rad, then our entire observable universe lies well-within this patch at the onset of 
inflation and is therefore highly homogeneous. 



4.3 Flatness and Homogeneity Problems in Ekpyrotic Cosmology 

Ekpyrotic theory solves the flatness and homogeneity problems in an analogous way, yet 
using drastically different dynamics. Instead of exponentially rapid expansion with nearly 
constant Hubble radius, the ekpyrotic phase consists of a slowly contracting universe with 
rapidly shrinking Hubble radius. The scale factor is essentially constant during this phase, 
as indicated in Fig. |4j Meanwhile, the Hubble radius shrinks by an overall exponential factor 
gA^ek^ with iVek defined by 

^^■^ = n n ]-\n[ . (4.6) 

This definition implies that e^'^'^ is just the ratio of the comoving Hubble radius at the onset 
and exit of the ekpyrotic phase. Thus by the end of the ekpyrotic phase the curvature 
component is suppressed by 

O(ck-cnd) / rr \ 2 

"fc _ / "ck-bcg-nck-bcg \ _ g-2Afek ('4 ^\ 

^(ek-beg) l^aek-endi^ek-end/ 

This is followed by a bounce and an expanding hot big bang phase. In the New Ekpyrotic 
scenario, this is achieved by a ghost condensation phase that violates the NEC, generates 
a non-singular bounce and eventually reheats the universe. As indicated in Fig. [4], how- 
ever, between tek-end and treheat the scale factor changes by at most a factor of order unity. 
Meanwhile, although the Hubble radius changes dramatically in the process — H vanishes 
momentarily at the bounce by definition — its magnitude at the beginning and end of the 
NEC- violating phase is nevertheless essentially the same, as shown with explicit bouncing 
solutions [IH]. Ignoring the details of this intervening phase, therefore, we can safely assume 
that a and \H\ both match continuously from the exit of the ekpyrotic phase to the onset of 
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the radiation-dominated expanding phase, so that 

^(ek-cnd) _ ^(reheat) _ ^^_g^ 

Therefore, just as in inflation, starting with ~ the onset of the ekpyrotic 

phase, one obtains an acceptably large radius of curvature by the present time provided 
^ck ^ ^rad- In this preclse sense the ekpyrotic scenario solves the flatness problem. 

The condition Nc-^ ^ A^rad implies that our entire observable universe is initially well 
within the proto-ekpyrotic patch. Therefore, the assumption that the initial Hubble radius 
is reasonably smooth guarantees a high degree of homogeneity for our observable patch. 

4.4 Observational Constraints on Model Parameters 

We have seen that inflation and ekpyrosis are equally successful at addressing the well-known 
problems of standard big bang cosmology. In particular, for the flatness problem, both allow 
for the natural choice Qk ~ initially. Of course, in absolute terms the initial radius 

of curvature is much larger in ekpyrosis than in inflation. As stressed earlier, however, 
this is a consequence of strikingly different assumptions about initial conditions in the two 
models: inflation proposes a hot beginning, corresponding to a microscopic proto-inflationary 
patch, whereas ekpyrotic theory prefers a cold and vacuous initial state, corresponding to a 
macroscopically large Hubble radius. 

To make this point quantitative, in this subsection we estimate the initial Hubble radius 
in the ekpyrotic scenario for a range of reheating temperatures. First, using A^ek ^ A^rad, we 



find from (4.6) that 



^0 

|-^ek— bcgl ^ ^ |-^ek— end| Ttt |-^ck— end| ; (^'9) 

reheat 



where in the last step we have substituted (4.2 ). As mentioned before and as shown in Fig. 4 



the net change in from the end of the ekpyrotic phase until reheating in the expanding 
phase amounts to a factor of order unity, and thus |-ffek-end| ~ -f^reheat- Moreover, from 
the Friedmann equation in a radiation-dominated universe, we have ifreheat ~ ^reheat/^Pi- 



Substituting this into (4.9), we obtain 



|i/ek-beg| <^^§^. (4.10) 

For Treheat = 10^^ GeV, this gives |ifek-beg| ^ 10^^ eV, corresponding to an initial Hubble ra- 
dius on the order of 1 m. For Treheat = 10^ GeV, we find |-ffek-beg| ^ 10"^'^ eV, corresponding 
to a Hubble radius of the order of 10^ km. 
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Note that the steepness of the V(0) factor in (3.6) is characterized by the size of the 
fast-roll parameter e defined in (3.7). In the more specific potential given in (3.2), where 
the V(0) factor and the quadratic term are the same as in the two exponential case, we find 
that the steepness of the exponential is controlled by the constant p = 2e. Restricting our 



discussion to (3.2 ), we note for future reference that the parameter p depends on the reheating 
temperature through the WMAP constraint on the amplitude of density perturbations — 



see (3.26) and (3.29). Neglecting the model-dependent prefactor P and again using p = 2e, 



we can solve (3.26) and (3.29) for p: 

p ~ 10^° 



reheat 



Mpi J 



10 



10 



reheat 



Mpi ) 



(4.11) 



For instance, Treheat = 10^^ GeV yields p ~ 10"^. Notice that the fast-roll condition p <^ 1 
combined with the WMAP amplitude puts an upper bound on the reheat temperature 
of Treheat ^ 10^^ — 10^^ GeV. This is no different than in inflation since, as mentioned 



earlier, (3.29) is identical to its inflationary counterpart. 



5 Initial Conditions and Taming the Instability 

In this section, we return to the tachyonic instability and study the associated initial condi- 
tions. We show that in order for the field trajectory to remain near x ~ = sufficiently 
long to generate the appropriate range of scale-invariant perturbations, x must be exponen- 
tially close to the top of the potential ridge at the onset of the ekpyrotic phase. Were this to 
be an initial condition, it would imply that ekpyrotic theory requires incredible fine-tuning 
to achieve a realistic cosmology. 

Happily, as we will show below, this condition can easily and naturally be satisfied without 
fine-tuning. This is accomplished by introducing a positive mass term for x, letting x couple 
to light fermionic degrees of freedom and allowing for a pre-ekpyrotic phase where this 
positive mass dominates. Hence, in this early phase, the system will predominately oscillate 
around x = 0- For a wide range of initial conditions, the energy in these initial x oscillations 
is efficiently transfered to the massless fermions, thereby bringing x exponentially close to 
= by the onset of the tachyonic phase. Meanwhile, the energy density in the radiation 
produced quickly becomes subdominant compared to that in 0, thanks to the ekpyrotic 
attractor mechanism. 

For ease of presentation, we will carry through our analysis not with the most general 



potential (3.6) but, rather, with the more restricted potential (3.2). We do this only for 
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convenience. Our results easily extend to the general case (3.6). Recall that an essential 
feature of (3.2) is that the steepness of the pure exponential potential in is determined by 
the parameter p <^ 1. 



5.1 Constraints on Initial Conditions 



To begin, let us carefully derive the constraints on the initial conditions associated with the 



instability in x during the ekpyrotic phase. From the perturbation analysis in Sec. 2.2, it 



is clear that the tachyonic instability is crucial in generating a Harrison-Zeldovich spectrum 
for 6x- This follows from the substitution of V^^ = 



into the fluctuation 

equation (|2.17|). Hence, we cannot remove this instability without spoiling scale invariance. 



tachyon 



The main constraint associated with this instability is the following: the initial value of x 
must be sufficiently close to Xt = so that the ekpyrotic phase driven by cf) lasts long enough 
to generate perturbations over the desired range of modes. 

A necessary condition for this to be the case is that the energy density in x be subdom- 
inant throughout the entire ekpyrotic phase. The total energy density during this phase is 
given by 

Ptotal = P^i + Px ' (5-1) 

where we have defined 



1 . 



1 



2 tachyon 



+ o{x' 



(5.2) 



with Vcff(0) and mtachyon given in (2.12) and (2.14), respectively. Hence, one must impose 
the constraint 

\Px\ ^ \pA (5.3) 
for all times during the ekpyrotic phase. In this case, we have 



Pij) ~ Ptotal 



(5.4) 



where the last equality is just the Friedmann equation. Using expression (2.16) and the fact 



that p <^ 1, it follows that one can ignore the Hubble term in the x equation of motion. 



From (2.11) and (2.14), this equation becomes 



whose growing-mode solution is 



X ^2 



X 



-1 



(5.5) 



(5.6) 



25 



It follows that the energy density in x is 



Px 



^ A ^tS'ChyonX 



(5.7) 



Returning to condition (5.2), we see that this constraint is strongest at tek-end marking 



the end of the ekpyrotic phase. At this time, it follows from (5.3), (5.4) and (5.7) that 



^^tachyon(^ck— cnd)^Xek— end ^ l-f^ck-endl M-pi 



where Axek-end is the displacement of the field x from the top of the ridge at the end of the 



ekpyrotic phase. Using (2.14) to substitute for mtachyon, we obtain 



AXck-cnd < P Mpi 



(5.9) 



Let us see what this entails for the x displacement at the onset of the ekpyrotic phase. 



denoted by Axet-beg- Noting from the scaling solution (2.11) and (5.6) that Ax ~ \H\, we 
find 



AXek-beg ^ 



H, 



ek— beg 



ek— end 



p M] 



PI 



p M] 



PI 



(5.10) 



where in the last step we have substituted (4.6). In other words, x must be exponentially 
close to Xt = at tek-beg in order for the ekpyrotic phase to last sufficiently long. 

We can also express this bound as a constraint on the initial energy density in x at the 
beginning of the ekpyrotic phase: 



(ek-beg) 

Px 

(ek-beg) 



'^tachyon (^ek-beg) AXgjj._,^gg 



< 



-"ek-beg^^-'Pl 
P '^tachyon(^ek-beg) 
-"ek-end 



ek-beg 



ek— end 



(5.11) 



Thus the initial x energy density has to be exponentially suppressed with respect to that of 
at the onset of the ekpyrotic phase. 



We conclude from (5.10) and (5.11) that the value of x at tek-beg has to be exponentially 
close to the top of the tachyonic ridge in order to generate the necessary e-foldings of scale- 
invariant perturbations. As an initial condition, this introduces a fine-tuning into ekpyrotic 
theory. However, as mentioned in the introduction to this section, the same result can be 
obtained naturally, without fine-tuning. We now show how this can be achieved. 
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5.2 Stabilization Mechanism 



In this subsection, we argue that the bounds (5.10) and (5.11) are naturally satisfied if the 



ekpyrotic phase, during which x is tachyonic, is preceeded by a phase in which x has positive 
mass squared and couples to light fermions. This can be achieved, for example, by adding 
a mass term m^(0)x^ to (3.2) where, in general, depends on 0. The potential then is 
given by 



V{<P,x) = -V,e 



1 + 



m\{4>)x^ 



The new mass parameter m^{(f)) is chosen so that effective mass term for 



(5.12) 



(5.13) 



is positive at early times but becomes tachyonic for sufficiently small values of 0. This 
can be the case even for constant m-^, since mtachyon grows in time as evolves down the 
exponential potential. However, this transition is facilitated by choosing m^{(f)) to be a 
decreasing function of 0. The end of the pre-ekpyrotic phase and the beginning of the 
ekpyrotic phase occur by definition at time tek-beg when x becomes tachyonic: 



^x(^ek— beg) ^ '^tachyon (^ek— beg) 



\H. 



ek— beg 



P 



(5.14) 



After this time the dynamics is dominated by the ekpyrotic potential (3.2). A typical po- 



tential of the form (5.12) is shown in Fig. 2 



In order to generate the desired number of e-folds of perturbations, must become 
subdominant at a sufficiently early time. Otherwise, the tachyonic phase will be triggered 
too late in the evolution of the universe, resulting in too short a range of scale-invariant 
perturbations. Assuming that the largest observable mode exits the horizon at tek-beg? a 



conservative assumption, we can use (5.14) to put an upper bound on the allowed value of 



at tek-beg- Expressing the right-hand side of (5.14) in terms of the reheating temperature 



by substituting for -ffek-beg and p through (4.10) and (4.11) respectively, we find that 

-10 ^0 f Mpi \^ ^ ^ _ -,^-40 f ^Pl 



xl''ek-beg. 



< 10- 



PI 



reheat 



Ml 



PI 



10 



reheat 



Ml 



PI 



(5.15) 



In the last step, we have used the fact that the microwave background temperature is Tq ~ 
2.7 K ^ 10-30Mpi. For example, T^eheat = 10^^ GeV gives m^(tek-beg) < 10"^ eV, whereas 
for Treheat = 10^ GcV the bound greatly weakens to ^^(tek-beg) ^ 10^ GeV. Finally, for 
Treheat ^ 10^ GcV, wc get a trivial bound of m^(tek-beg) ^ 0.1 Mpi. 
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The pre-ekpyrotic phase occurs during the time interval ti < t < tck-bog? where tj is 
some as yet unspecified initial time. By definition, during this phase x has positive mass 
squared given, for most of the epoch, by m^g ^ m^. Given some initial displacement A^j, x 
will undergo oscillations around the minimum of the potential at x = 0. Since the universe 
is contracting, we cannot rely on Hubble friction to damp out these oscillations. On the 
contrary, they will instead get amplified by the usual blueshift effect. However, if we assume 
that X is coupled to light fermionic degrees of freedom, then its coherent oscillations will 
decay into these particles, thereby stabilizing x- Let \l/ be a typical light fermion and assume 
that it couples to x ^-s X^^ ■ Then the rate of decay is given as usual by [ST] 

r-m^. (5.16) 

If r ^ the decay is rapid on a Hubble time and the x evolution is overdamped. That is, 
the field settles at its minimum well before the blueshift effect can amplify the oscillations. 



This is easily satisfied here if we choose at tek-beg to saturate the bound (5.14) since, in 
this case, we have 

r«x(^ek-bcg) ~ ^ |if^k-beg| Z\Hi\, (5.17) 

where Hi is the initial Hubble parameter. Since is either constant or increases at earlier 
times, it follows that F ^ \H\ at any time during the pre-ekpyrotic phase. 

Let us now calculate the initial conditions at the beginning of the pre-ekpyrotic phase that 



will produce the exponentially suppressed displacement (5.10) and ratio of densities (5.11) 
at tek-beg- As in the ekpyrotic phase, we will assume for concreteness that 

so that the evolution is dominated by during the entire pre-ekpyrotic epoch. It will soon 
become obvious that our stabilization mechanism is sufficiently effective to allow for larger 
initial x component. Of course the mechanism eventually fails if is initially orders of 
magnitude larger than p^, but these are clearly fine-tuned and unnatural initial conditions. 
To simplify the calculation, we henceforth take to be independent of 0. Moreover, 



we assume that saturates (5.17), 

~ "^tachyon(ick-beg) (5.19) 

so that r ^ \H\. In this regime, we can ignore any blueshift effect. It follows that from the 
earlier discussion that 



1 
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m^Ml, , IpJ ^ -mix'- (5.20) 



Putting these densities into (5.18 



) and evaluating the inequahty at gives 



m^{ti)Ax^ < WMpi , (5.21) 
where Axi is the displacement of field x the beginning of the pre-ekpyrotic phase. Us- 



ing (5.17) and (5.19), we find that the initial displacement must satisfy 



Ax^<{J^^]pMpl. (5.22) 



ek— beg 

This consistency condition ensures that x is subdominant throughout the pre-ekpyrotic 
phase. 

To determine the allowed values of Axi, one must evaluate the ratio ifj/ifek-beg- This is 
most easily done by considering the ratio of the energy densities. From the above discussion 
it follows that the energy density in x decays as oc e~'"* ~ e~^^^. Thus, starting with 
some energy density p^^ ^ii by the onset of the ekpyrotic phase this has decayed to a value 

^^ek-beg) gj^gj^ 

(ek-bcg) 

~ g-nix(*ck-bog-ti) r:^ Q-m^i-U) ^ g-V2|//ck-bcg|(-ti)/P 23) 

(i) ^ ~ ~ ) V • / 

Px 

where in the second step we have assumed, for simplicity, that the pre-ekpyrotic phase is 
long compared to tek-beg; that is, —U ^ — tek-beg- Moreover we have substituted ~ 



"^tachyon(4k-beg) ~ | -f^ek- beg | /P frOm (2.14). 



It is convenient to rewrite this in terms of the initial Hubble parameter Hi. Since 
Px"* ~ P^6 by assumption, the cosmological evolution is dominated by throughout the 



pre-ekpyrotic phase. Therefore, we can substitute the relation Hi = p/ti given in (2.11) 



and thus 



(ek-beg) 

^ e-V2^^ok-bcg/J^« ^ (5 24) 

Px 

(ek— beg) (i) 

P]^ ^ g-V^gck-bcg/g» ^ g-V2^^ok-bog/J^» _ (5.25) 

(ek-beg) (ek— beg) ' \ ' / 

P(j, Pcf) 

The last inequality follows from p^^^ ^'^^^ > p^^ (due to blueshift between U and tek-beg) as 
well as p^^ > p^x (by assumption). Comparing this against the bound given in (5.11), we 
see that it is satisfied if 

^^>v^iVek. (5.26) 
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Putting this into (5.22), we get an expression for the displacement at the beginning of the 

1 



pre-ekpyrotic phase given by 



PI - 



(5.27) 



To get a feel for what this expression entails, recall from Section 4 that p and Ne]^ depend 
on the reheat temperature. For example, for Tj-eheat = 10^^ GeV, we found that p ~ 10^^ 
and Ne]^ > 60. It follows that for this reheat temperature the displacement at the beginning 
of the pre-ekpyrotic phase is given by 



Ax. < 10 



(5.28) 



representing a wide range of natural initial conditions. 

It is also useful to consider the ratio of the displacement of x ti to the exponentially 
fine-tuned value of x required at tek-beg- Using (5.10) and (5.27), we find that 

. (5.29) 

AXek-beg V^A^ek 

Clearly our mechanism naturally removes the exponential fine-tuning problem, as claimed. 
The remaining issue is to understand the physical implication of choosing the Hubble 



parameters at tek-beg and to satisfy the ratio (5.26). Using the fact that during the 

: follows 

>v^iVek, (5.30) 



^-dominated pre-ekpyrotic phase H = p/t, it follows from (5.26) that 

U 



^ek— beg 

a natural and physically reasonable assumption. For example, for the reheat temperature 
^reheat = 10^^ GcV this ratio becomes 

U 



>10\ (5.31) 

lek— beg 

We conclude that in order for x to be close enough to the top of the tachyonic ridge 
by the onset of the ekpyrotic phase to ensure a sufficient number of e-foldings, it suffices to 



start out at ti with the large, natural value of x given by (5.27). This will be the case as 
long as the pre-ekpyrotic phase begins sufficiently before the onset of ekpyrosis, the precise 



relationship being specified in (5.30). 



6 Graceful Exit and Bouncing 

Thus far, we have discussed 1) the pre-ekpyrotic phase in which x gets stabilized near x* = 
and 2) the ekpyrotic phase in which the entropy perturbation acquires a nearly scale-invariant 
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spectrum. In this section, we turn to the final two phases required for a reahstic cosmology, 
that is, 3) the exit from the ekpyrotic phase and 4) the bounce from contraction to expansion 
enabled by a NEC-violating ghost condensate. We begin with a general discussion of these 



last two phases, in Sees. |6.1| and |6.2| respectively, before turning to a specific example in 
Sec.m 



6.1 Converting Entropy into Curvature Perturbations 



The ekpyrotic phase described by the scaling solution (2.11) must eventually come to an end 
if the universe is to undergo a smooth bounce and enter the hot big bang expanding phase. 
Before turning to the bounce, however, we first discuss a key element of New Ekpyrotic 
Cosmology which also occurs during the exit from ekpyrosis, namely the conversion of the 
scale-invariant entropy spectrum into the adiabatic mode. Even before the bounce, the 
exit from the ekpyrotic phase can be used to imprint the scale-invariant entropy spectrum 
onto the adiabatic mode, described by (, the curvature perturbation on uniform-density 
hypersurfaces. 

In general, it is convenient to consider new field variables, a and s, defined with respect 
to the field trajectory. As shown in Fig. [5} the adiabatic field velocity, &, is proportional 
to the vector tangent to the curve, while the entropy field velocity, s, is proportional to the 
normal vector. In other words, their decomposition in terms of (f) and x is given by 

& = cos 9 (p + sm9 X 

s = — sm9 (p + cos9 X , (6.1) 



where 



tan^ = ^. (6.2) 



During the ekpyrotic phase, the background field motion is along 0, and thus 9 = 0. 



Evidently, it follows from (6.1 ) that in this phase a and s coincide with (j) and respectively. 



Once we exit the ekpyrotic phase, this is of course no longer true. In general, the (0, x) 



equations of motion given by (2.10) imply the following evolution equation for the adiabatic 
field 

a + 3H& = -V„ , (6.3) 



where Vcr = cos9V^ + sin6'V^^. This result will come in handy later on in the discussion. 
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X 




Figure 5: Adiabatic (a) and entropy (s) directions at a given point along the field trajectory 
in (0, x) space. 



The relation between the adiabatic and entropy fluctuations to 50 and 5x follows imme- 



diately from (6.1 ): 



Sa = cos 9 S(f) + sin 9 6x 
6s = — sin 9 6(f) + cos 9 6x ■ 



(6.4) 



Using the standard energy and momentum constraints, the entropy perturbation satisfies |52j 



6sk + 3H6sk + {^ + Vss + S9'' ] 6sk 



2 G k'^ 



(6.5) 



cos^ 9 Vy.y.+sm'^ 9 



where is the curvature perturbation in Newtonian gauge, while Vg 
is the curvature of the potential orthogonal to the field trajectory. 

Since 9 = during the ekpyrotic phase, as mentioned earlier, we have 6s = 6x, so that 
fluctuations in x indeed correspond to entropy perturbations. As a quick consistency check, 
since ^ = and Vss = Kxx' ^^^^ (6-5) indeed reduces to (3.10). Note, however. 



that the identification of 6x with the entropy perturbation only holds during the ekpyrotic 
phase; that is, as long as the field trajectory is along the direction. We will shortly 
introduce a feature in the potential that causes a sharp turn in the field trajectory at the 
end of the ekpyrotic phase, thereby imprinting a scale-invariant spectrum on (, which is 
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proportional to the adiabatic fluctuation 5a. Subsequently the classical field evolution is 
partially along x, and, evidently, fluctuations in the latter are no longer identifled with the 
entropy perturbation. 

Let us therefore turn to the conversion of the entropy fluctuation to the curvature per- 
turbation spectrum. The latter is sourced by the entropy perturbation as follows 



C,k = ~'^k + —e5sk. (6.6) 

The flrst term proportional to fc^ is seemingly negligible at long wavelengths. As pointed out 
in [33], however, this is not necessarily the case during the exit from the ekpyrotic phase. 
For our purposes we will neglect this term altogether for simplicity. 

Clearly, to endow with a scale-invariant contribution from 5s, the fleld trajectory must 
undergo some turn at the end of ekpyrosis, which we assume is sharp and therefore rapid 
on a Hubble time. In [18], this was achieved by introducing a minimum followed by a sharp 
rise in the respective potentials for 0i and 02. In the simplifled 0, x framework studied 
here, this would correspond to adding new terms in the potential that depend on speciflc 
linear combinations of and x- A much simpler option suggests itself, however, namely to 
introduce some term which pushes x away from = at the end of the ekpyrotic phase. 



To be precise, we want to modify potential (3.2) so that = is only an approximate 
solution to the equations of motion. One can, for example, add a term to the potential 
which is approximately linear in x near x = and therefore prevents the latter from being 
an extremum of the potential. This linear term must be a function of so that it is negligible 
during the ekpyrotic phase but eventually becomes important after a sufficient number of 
e-foldings is achieved. At this point it tips x away from the tachyonic ridge marking the end 
of the ekpyrotic phase. This is sketched in Fig. |3j 



To illustrate this, let us describe the ekpyrotic phase by the simple potential (5.12) and 
add a correction of the form 

0/Mpi ( ^ , ]_ _ ^2 



= t/oe-V ^7-/-^- 1 + — ^(X - Xo)^ . (6.7) 



Clearly, this has non-vanishing slope at x = 0. For p' <^ the exponential prefactor makes 
this correction vanishingly small during the ekpyrotic phase, and therefore Xt ~ is still a 



valid solution then. Once becomes sufficiently small, however, the correction (6.7) becomes 
important and pushes x away from the x* = tachyonic point. This example will be studied 
in detail in Sec. [7l 
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Coming back to general considerations, we assume that the driving force is sufficiently 
strong to generate a fast-roll of x away from the tachyonic point, corresponding to a sharp 
turn in the field trajectory. To make this more precise, we require that the driving term V^^ 



in the equation of motion for x — see (2.10) — is comparable to the Hubble damping term: 



lKx(x = 0)| >Mpiif2. (6.8) 
Equivalently, this condition amounts to requiring the slow-roll parameter = y^/if^M|[ to 



be comparable or greater than unity. For example, for the specific correction given in (6.7) 



this condition is satisfied for sufficiently small values of (j). We henceforth denote by 0ek-end 



the value of when (6.8) is fulfilled. This notation is consistent with Fig. |4| since the rolling 
of X away from the tachyonic ridge marks the end of the ekpyrotic phase for all practical 
purposes. In particular, fluctuation modes produced thereafter are no longer scale- invariant. 

In the limit of fast-roll, we can treat the turn in the field trajectory as instantaneous on 
a Hubble time. This amounts to approximating 9 with a delta function and treating H as 
constant. Hence 

Cfc ~ ^l^Z^A^ 5{t- tek-end) , (6.9) 

a 

where AO is the total change in angle during the exit. Of course, smoother exit dynam- 
ics could do equally well, but this rapid-exit approximation makes it possible to proceed 
analytically. 

As argued in [18], a is constant during the instantaneous exit from the ekpyrotic phase. 



This follows by inspection of the equation of motion for a given in (6.3). If a is not constant 
during the transition, then the '& term will generate a delta-function contribution. But this 
cannot be compensated by any other term in the equation. The 3if(T evidently cannot do the 
trick, since H is constant during the transition by assumption. Similarly, a sudden change in 
field direction can at most yield a jump in V^rri but not a delta-function discontinuity. Thus 
a is continuous, and we can substitute its value at the end of the ekpyrotic phase. Since 



a = (f) during this phase, we read off from (2.11 ) 



a = = ^ — Mpi , (6.10) 

t-ek— end V £ 

where we have used p = 2e. 

Similarly Ss changes by at most a factor of unity during the exit, and we refer the reader 
to [IH] for a detailed argument. Thus we approximate 6s to match continuously as well, for 
simplicity, and can substitute for 6s its value at the exit from the scaling solution. Since 

34 



Ss 



Sx up to the exit, this can be read off from (2.11 ) and (2.19), again using p = 2e: 



\H, 



ck"Cnd 



23/2( 



(6.11: 



Substituting the above expressions for & and 6s, we can integrate (6.9) and obtain 



'ek— end | 



A^ 



reheat 



(6.12) 



where we have used \H, 



ek— end 



/2eMpi V2eMpi 
-f^reheat indicated in Fig. liland discussed earher. Compar- 



ison with (3.29) shows that the model-dependent factor P is identified with A9 in our sharp 



turn approximation, thereby confirming the discussion below (3.29). It is desirable to have 
AO ~ C'(l), for otherwise one needs even smaller values of e at fixed reheating temperature 



to maintain the WMAP normalization (3.26). In the next subsection we will discuss the 
conditions for which AO ~ 0{1). 

To summarize, the end of the ekpyrotic phase is set by a term in the potential which 
pushes X away from the Xt = ^ tachyonic point. This term is assumed irrelevant during 
the ekpyrotic phase, so that x ~ is a good approximate solution, but eventually becomes 
important and triggers the end of ekpyrosis. This sudden motion in x corresponds to a 
turn in the field trajectory, which in turn imprints the scale-invariant entropy perturbation 
spectrum onto the curvature perturbation. 

Instead of dialing the end of the ekpyrotic phase with additional terms in the potential, 
one could alternatively exploit the tachyonic instabihty in x, as studied in [321 133]. However, 



as the discussion of Sec. |5.1| most emphatically underscores, to generate the desired number of 
e-foldings in this case requires an exponential fine-tuning in the initial x displacement — see. 



for example, (5.10). If the initial displacement is larger, then the ekpyrotic phase will be too 
short, resulting in an unacceptably narrow spectral range of scale- invariant perturbations. 
If it is smaller, on the other hand, then the ekpyrotic phase will last for too long. That 
is, by the time the entropy perturbation is converted onto (, the Hubble parameter will be 
too large in magnitude to yield an acceptable perturbation amplitude. In other words, the 
WMAP constraint on the amplitude of perturbations is tied in this latter case to a tuning 
of initial conditions. 

On the other hand, adding correction terms to the potential, as we do here, leads to robust 
predictions for the perturbation spectrum for a broad range of initial conditions. Moreover, 



we know that corrections to (5.12) are necessary for independent reasons. For instance, 
the NEC-violating ghost condensation phase requires a nearly fiat and positive potential. 
Thus the negative, steep exponential potential of the ekpyrotic phase must eventually have 
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a minimum and rise to positive values. What we are advocating is that these corrections will 
generically push x away from x = 0, thereby endowing Q with a scale-invariant spectrum. 



We will illustrate this in detail in Sec. [7] with the correction (6.7), showing how it leads 
to a scale-invariant curvature perturbation and how it can be matched smoothly to the 
NEC- violating ghost condensate phase. 

6.2 A Bouncing Scenario 

In this subsection, we complete the scenario by discussing the dynamics of a non-singular 
bounce along the lines of [IH]. This is achieved by supplementing the ekpyrotic phase with 
a phase in which the energy density is dominated by that of a ghost condensate 120]. The 
ghost condensate can violate the null energy condition [23], a necessary ingredient to produce 
a bounce, without introducing ghost-like instabilities. The analysis closely parallels that 
of [IH], and thus we will be brief, referring the reader to our earlier paper for details. A 
key difference, however, is that here it suffices to have only one field entering the ghost 
condensate phase, namely the ekpyrotic field 0; meanwhile x is a standard scalar field with 
canonical kinetic term throughout. 

i) Rolling of x towards minimum 

Before entering the ghost condensate phase, we assume that x g^ts stabilized. As de- 
scribed in the previous subsection, the production of scale-invariant modes comes to a halt 
once corrections to the ekpyrotic potential drive x away from the tachyonic point. Such 
corrections are assumed negligible at early times but become important for smaller than 
some critical value 0ek-end, where the subscript indicates that this marks the end of the 
ekpyrotic phase for all practical purposes. 

Subsequently we assume that x rolls towards a stable minimum, denoted by Xmm- For 



instance, our fiducial correction term given in (6.7) generates a minimum at Xmin ~ Xo- Note 



that the fermionic degrees of freedom which stabilized the field in the pre-ekpyrotic phase 
now become heavy as x acquires a non-zero expectation value. Thus they become irrelevant 
and cannot be relied upon to stabilize x this post-ekpyrotic phase. Instead x rolls down 
and underg oes undamped harmonic oscillations around Xmm- 

The kinetic energy acquired by x as it rolls away from the ridge is of course just given 
by the difference in potential energy between x = and x = Xmin- 

Px(^ek-cnd) = 2 ^^(^ck-cnd) = V (^ck-cnd, 0) — V (^ek-end, Xmin) • (6.13) 
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If we assume for simplicity that x is massive compared to Hubble, 



(Xmin) » , (6.14) 

then the field undergoes many oscillations in a Hubble time. Consequently, its energy density, 
averaged over many oscillations, blueshifts as dust: 

p^-a-\ (6.15) 

Although this grows in time as the universe contracts, nevertheless quickly becomes 
negligible compared to the energy density in 0, which keeps rolling down its steep quasi- 
exponential potential. In other words, thanks to the ekpyrotic attractor mechanism, we can 
tolerate ^ p^f, at tek-end since (f) will dominate again within a few e-folds of contraction. 
These considerations impact on the allowed range of A^, the change in angle in the field 



trajectory as we exit the ekpyrotic phase. As shown in (6.12), this parameter enters in the 
amplitude of density perturbations. Here we argue that a large change in angle is indeed 
possible. 



From the definition of 6* in (6.2), we see that A6' ~ 0{1) is achieved provided 



X (^ck-ond) ^ (i^ck-end) . (6.16) 



Equivalently, since p^j,- (tek-end) ~ X (tek-end) /2 from ( |6.13[ ) whereas p^K. p (j)^/2 from ( |2.1l[ ), 
we need 

Px'> P<P at t = tek-end . (6.17) 

The discussion of the previous paragraph shows that such a large ratio of densities is allowed 
by the ekpyrotic attractor mechanism. In other words, x can dominate momentarily as it 
rolls off the ridge and starts oscillating, thereby generating a significant AO. Very soon, 
however, will take over as it keeps rolling down its steep negative exponential potential. 

Once X is subdominant, the dynamics effectively reduce to that of a single field 0, and 
the story is virtually identical to single-field ekpyrosis with NEC-violating ghost condensate 
and non-singular bounce. The effective single-field potential is shown in Fig. |6} Region a) 
describes a steep negative exponential potential corresponding to the ekpyrotic phase. A 
necessary condition to violate the NEC is that the potential becomes positive. To see this, 
recall that the onset of NEC- violation occurs when H = 0. And since Mp^H = —0^/2, 
this also corresponds to 0^ = 0. But then the Friedmann equation, 3H^Mp^ = 0^/2 + V, 
immediately tells us that V > 0. Thus the field must reach a minimum in V and rise 
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to positive values (region b)). Such a minimum and sharp rise can be achieved with our 
ekpyrotic potential (5.12) plus the correction (6.7). The region c) corresponds to the ghost 
condensation phase where the potential is approximately flat. 




ii) Ghost condensate phase and NEC violation 

Our next consideration is the epoch of NEC violation, triggered by entering a ghost 
condensate phase. This is achieved by invoking higher-derivative corrections to the kinetic 
term for (p, which takes the form 

C = v/=^M^P(X), (6.18) 

where 

X.^^^ (6.19) 

is dimensionless. Here m and M are some mass scales to be determined by the fundamental 
theory. As successful merger of ekpyrosis and ghost condensation, as we will see, requires a 
large hierarchy: M ^ m. 

The ghost condensate relies on the kinetic function P{X) having a minimum at some 
finite X, taken to be 1/2 without loss of generality, corresponding to = — m^. In the 
absence of a potential, X = 1/2 is an exact solution to the equations of motion, even in a 
cosmological background. 
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What about the global form for P? During the ekpyrotic phase, is assumed to have 
approximately standard kinetic term, corresponding to P{X) ^ X. In order for fluctuations 
to have positive norm, this quasi-linear part must lie to the right of the minimum, i.e. for 
X > 1/2 [18]. If the quasi-linear part lies at X < 1/2, on the other hand, the field evolution 
must go through a region with Px < 0, corresponding to ghost-like fluctuations [IH]. The 
desired form of P{X) is sketched in Fig. [7j 

P(X) 




Figure 7: Kinetic function for (j). 

In the process of climbing up the sharp rise in the potential, the kinetic energy in 
decreases dramatically, and X is brought to the vicinity of the minimum of P{X) by the 
time reaches the plateau. This marks the onset of the ghost condensate phase which for 
convenience we will set at t = 0. Then, near the ghost condensate point we can expand 
the field as = —mH + n, leading to the following expressions for the energy density and 
pressure: 

p = M\2PxX-P) + V^-^^ + V 

V = M^P{X)-V^ -V, (6.20) 

where K = P^xx- See [18] for details. (Notice that we have taken P(l/2) = without loss 
of generality since an overall shift in P is degenerate with a shift in V^.) 

The culprit in violating the NEC is the term linear in tt, since its contribution to the 
energy density does not have a definite sign. Put another way, the H equation is given by 

M^,if = 4(p + P)-^, (6.21) 
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which clearly can take either sign depending on the sign of tt. 



in) Constraints for successful bounce 

The successful merging of ekpyrosis with ghost condensation imposes some constraints 
on the ekpyrotic potential. 

Firstly, by assumption we must have X > 1/2 throughout the ekpyrotic phase, that 
is, 0^ ^ m'^. Because the kinetic energy grows during the ekpyrotic phase, evidently this 
condition is most stringent at the onset of this phase: 0^(tek-beg) ^ And since 0^/2 ^ 
—V during the ekpyrotic phase, this immediately implies: |\^(0ek-beg)| ^ w^"^- For our 
purposes we find it convenient to express this as a constraint on the value of the potential at 



the minimum, denoted by Vmin = ^(0min, Xmin)- From the scaling solution (2.11) and (4.6), 
we have |V^(0ek-beg)| ~ e~^^'='^|Vniin|- Therefore we must require 

iKiinI > e^^^m^ (6.22) 

Secondly, by assumption (f) must lie in the vicinity of the ghost condensate point by the 
time it reaches the plateau, so that we can approximate ~ —m?t + tt. As shown in pj] , 
this gives an upper bound on Knin, the value of the potential at the minimum: 

Klin < • 6.23 

P 

To summarize, the ekpyrotic field is successfully morphed into an NEC-violating ghost 
condensate provided |Knin| falls within the range 

e^^-m^«|Knin|«— . (6.24) 

P 

Evidently this is satisfied provided M is exponentially larger than m. For instance, for a 
reheating temperature of 10^^ GeV corresponding to A'ek ^ 60 and p = 10^^, we can take 
M > 10^6 GeV and m > 10^ GeV. 

So far everything we have said is a review of aspects of ghost condensate/ekpyrosis merger 
first discussed in p.8j. A key difference with the two-field scenario of [IH], however, is that 
here we shall use only one field, namely 0, as ghost condensate, while x is stabilized. While 
this greatly simplifies the picture, there is an additional constraint to consider, having to do 
with the energy density in the x-oscillations, denoted earlier by p^. This is a priori non- 



trivial since the latter blueshifts as as the universe contracts — see (6.15) — , whereas 
the ghost energy density decreases since it violates the NEC. Nevertheless, this is not a big 
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concern for two reasons. As mentioned earlier, a few more e-folds of rolling down its 
steep exponential potential after the end of the ekpyrotic phase makes negligibly small 
compared to p^f, by the onset of the NEC-violating era. Moreover, from the moment enters 
the ghost condensate phase until reheating, the scale factor a{t) changes by at most a factor 
of order unity. Thus p^ does not blueshift appreciably throughout the NEC- violating epoch. 

It remains to show that a{t) satisfies this property. Clearly this is the case from the end 
of ekpyrosis until reaches the plateau, since everything happens almost instantaneously 
on a Hubble time, by assumption. That this is also satisfied throughout the subsequent 
NEC- violating phase follows from requiring that the bounce occurs within at most a Hubble 
time or so, which, as argued in p3j, ensures that the Jeans-like instabilities intrinsic to ghost 
condensation remain under control during the bounce. A numerical solution obtained in jl8] 
and plotted in Figs. [8] and [9] shows that the scale factor changes by approximately a factor 
of 2 over the course of the NEC-violating phase. See the captions for details. Consequently, 
p^ blueshifts by a factor of 8 during this phase. 




Figure 8: Evolution of the scale factor a{t) during the bounce generated by the ghost con- 
densate. The NEC-violating phase begins by definition when H = 0, which as shown in 
Fig. [9] occurs around t ^ 1.2 in these units. Meanwhile the bounce occurs by definition when 
d = 0, which occurs around t ^ 2.4. The point is that in this time interval the scale factor is 
seen to change by a factor of 2, from approximately 0.4 to 0.2. Note that we have set t = 
to be the beginning of the ghost condensation phase. 

Meanwhile, it is easily seen that the ghost condensate energy density remains essentially 
constant during this phase. To estimate its value, we assume for simplicity that the scaling 



solution (2.11 ) holds all the way to the onset of the ghost condensate phase. In this case we 
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Figure 9: Evolution of the Hubble parameter H{t) during the bounce generated by the ghost 
condensate. The epoch of NEC violation starts when H = 0, which from the plot is seen to 
occur at t ^ 1.2. The bounce {H = 0) occurs at t ^ 2.4. Note that we have set t = to be 
the beginning of the ghost condensation phase. 



have 



Pghost ~ P I Vyii 



(6.25) 



Therefore, ignoring the effect of the ekpyrotic attractor mechanism between tek-end until tc, 
a sufficient condition to ensure the subdominance of during the bounce is 



ck— end / 



V {4>ck-end, 0) -V {(pck-end, Xmin) ^ 8p \V^ 



(6.26) 



where the factor of 8 comes from the blueshift of py. through the bounce, as described above. 
In other words, the potential drop experienced by x it rolls towards its stable minimum 
must be small compared to V^in- 

6.3 Summary 

Let us summarize the sequence of events from the end of the ekpyrotic phase through the 
non-singular bounce and onto the radiation-dominated hot big bang phase: 



The entropy field x is eventually driven away from the tachyonic ridge, thereby bringing 
to a halt the production of scale-invariant density perturbations. See Fig. |3| This is 



achieved by a correction to the ekpyrotic potential (5.12), for instance of the form (6.7) 



which has non-zero slope at x = 0. For simplicity, we assume that the departure from 
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X ~ occurs rapidly on a Hubble time, which requires a fast-roll condition to be 
satisfied 

|^x(x = 0)| >//^Mpi. (6.27) 
We denote by 0ek-end the value of at which \ Vy^\ starts being comparable to H'^Mpi. 

Subsequently x rolls to a stable minimum denoted by Xmin and undergoes oscillations. 
As X acquires a non-zero expectation value, the light fermions of the pre-ekpyrotic 
phase become heavy and cannot be relied upon to stabilize x- We also impose that x 
is heavy at its minimum, 

T<Xx(Xmin)»//', (6.28) 

SO that its oscillations arc rapid on a Hubble time. The energy density acquired by 
X, given by P;^(tck-cnd) = V (0cnd,O) - V (0end,Xmm), blucshlfts effectively as dust and 
quickly becomes subdominant to (f). 

The evolution then reduces to single-field ekpyrosis. The end of the ekpyrotic phase 
is triggered by reaching a minimum in the potential, which subsequently rises to 
positive values where it becomes fiat. We denote by Kiin the value of the potential at 
the minimum, and by ^min the corresponding field value. Since stabilization of x is 
assumed to occur before hits the sharp rise, for consistency we must have 

0min ^ 0ek-end ■ (6.29) 

As climbs up the sharp rise of the potential, its kinetic energy decreases substantially, 
which brings X to the vicinity of the ghost condensate point. This marks the beginning 
of the NEC- violating phase. A successful synergy of ekpyrotic and ghost condensate 
cosmology requires V min to fall within the allowed window 

e'^^-^m' < V^i^ < , 6.30 

P 

where m and M parametrize the location of the minimum and curvature of the kinetic 
function P{X), respectively. 

In order to have a bounce, the energy density in the ghost condensate, given by pT^in, 
must dominate over the energy density in x oscillations, p^. Since the latter blueshifts 
by a factor of 8 or so during the bounce, this leads to the constraint: 

V (0ek-end, 0) - 1/ (0ek-end, Xmin) ^ 8p |Kiin| • (6.31) 
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These five equations are the main constraints to have a successful exit from ekpyrotic 
phase followed by a non-singular bounce. In the next section, we illustrate how these con- 
ditions can be fulfilled by studying a specific potential consisting of the ekpyrotic poten- 



tial (5.12) with the fiducial correction term (6.7). 



7 A Specific Exit Model 

We now illustrate the constraints described in the previous section with an explicit choice of 



potential. Our fiducial potential consists of two terms: i) the ekpyrotic potential (3.2) which 



we assume has exact exponential form in 0; ii) the correction term (6.7) which generates a 
minimum in (p and forces x away from the tachyonic ridge. Thus the full potential is given 
by 

mx) = Kk(0,x) + AV(0,x), (7.1) 

with 



V..{<P,X) = -^oe-V.*/*^- ( 1 + + . . . 



'PI 

AV{^,X) = f/oe-V7<^/*^- ( 1 + _L_(;^ - + ...), (7.2) 



where Vq and Uq are both positive, and, as earlier, ellipses denote higher-order terms in 
X which are irrelevant for our purposes. Note that we neglect the mass term for x added 



in (5.12). While this term is crucial in stabilizing x early on, it becomes subdominant by 
the onset of the ekpyrotic phase and remains subdominant forever after. The first term 
is the most constrained part of the potential since it pertains to the generation of density 
perturbations, whereas there is considerable freedom in the correction term AV. We focus 
here on this simple form for the latter given its manifest resemblance to the ekpyrotic piece. 

The exponents satisfy p' <^ p, which guarantees that AV is negligible during the ekpyrotic 
phase where (p is large. As (p decreases, however, eventually this term becomes relevant and 
pushes X away from the tachyonic point, triggering the end of ekpyrosis. 



7.1 Global Minimum for V 

This potential has a global minimum. Under some mild assumptions we can derive explicit 
analytical expressions for the corresponding field values (prnin and Xmin- First note that the 
condition dV/dcp = can be made nearly independent of x if we assume that the O(x^) 
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terms are small compared to unity, at least in the relevant range < x ^ Xo- (We will check 
a posteriori that Xmm ~ Xo so that this is indeed the range of interest.) Since p' <^ p it is 
sufficient to impose 

<V^. (7.3) 



Mpi 

And since p' ^ p ^ 1 this says that the minimum along x from the AV term lies at small 
values compared to Mp\. In other words, x moves a small distance in field space as it rolls 
away from the tachyonic point towards the minimum. We stress that this condition is not 
necessary but is only imposed to simplify the expression for 0min- We find 

. Mp, ,J^^^), ,7.4) 



Next we can solve dV/dx = for Xmin by substituting the above expression for (p^i^ 
everywhere: 

p Uo \ P 

It follows that AV^ plays the dominant role in determining Xmin- 



7.2 Exit from Ekpyrosis: x Starts to Roll 

The next step is to find when x starts to roll away from the tachyonic point at x ~ towards 
the minimum at Xmin ~ Xo- As described in Sec. 6.3, this occurs at a value 0ek-end at which 

\Vjx = 0)\r^Mp,H\ 



From (7.2) the left-hand side is given by 

vAx = 0) = - 



Une 



r0ok-cnd/Mpi 2 Xo 



p'Mpi 

Meanwhile, since the ekpyrotic phase is still going on until (p = 



scaling solution (|2.11|) to determine the right-hand side: 

pVoe 



-"ok-end^^^Pl 



Combining the above two expressions allows us to solve for (pcnd- 



M, 



PI 



^ek— end 



/p^ 



In 



2 ^0 Xo 



(7.6) 
(7.7) 

^ek-end, WC CaU USC the 

(7.8) 
(7.9) 



V2 y/p-VF'" \PP' Vo Mpi ^ 
To summarize, for > 0end, a good approximate solution for x is just x ~ 0. This is 
the ekpyrotic phase, characterized by rolling in the direction while x sits idle at the top of 
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the tachyonic ridge. Once reaches 4>end, however, we have \V^\ ~ H'^Mpi and x is pushed 
away from its unstable point. Subsequently we have \Vx\ ~ H^Mp\, thereby satisfying the 



first condition (6.27) of Sec. 6.3 



7.3 Mass of x about Minimum 

Once X reaches Xmin ~ Xo, it undergoes oscillations about this stable point. The second 
condition (6.28) ensures that the mass of x around Xmin ~ Xo is large compared to Hubble, 



so that the field undergoes many oscillations in a Hubble time. For simplicity we assume that 
(f) is nearly still as x rolls to the minimum, so that all relevant quantities can be evaluated 

&t 0ek-end- 



Using (7.9) we obtain 



,xx 



mm? Xinmj 



Mi; 



7</'ok-cnd M^Pl 



Ml 



PI 



(7.10) 

Xo ' 

Evidently this must be positive if we want Xmin to be a stable point by the time we reach 
0ek-end- The couditiou T^xxl'i^min) > amouuts to 

^<-- (7-11) 

Mpi 2 ^ ' 



Much like (7.3), this forces Xo to be small in Planck units. 



Using (7.8) we obtain 



.XX 



■Tniii) Xminy 



PI 



-"ck-cnd 



Xo 



(7.12) 



Barring some fine-tuning between these two terms, this ratio is generically much bigger than 
unity. It then follows that x acquires a large mass at the minimum and oscillates rapidly. 



7.4 Exit Happens Before Reaches Minimum 

Let us turn to the third condition, which requires that (p^i^i ^ <^ek-end- This ensures con- 
sistency of the sequence of events assumed here, namely that x starts to roll towards its 
minimum, thereby imprinting the entropy perturbation onto (, before reaches the steep 
rise in the potential. Again this assumption is not necessary — nothing prevents us from 
converting the entropy mode into the curvature perturbation after the has reached its 
minimum. However an estimate of ( would probably require numerical analysis in this case. 

in ^ 0end is Satisfied if 



In any case, comparison of (|7.9|) with (|7.4|) shows that 

'p' . 



Xo 



Ml 



PI 



(7.13) 
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Note that this condition is consistent with (7.3) and (7.11) since p' ^ p ^ 1. 



7.5 Energy Density in x Oscillations 



The last condition (6.31) puts an upper bound on the energy density in x it rolls to the 



minimum and starts oscillating. As argued in Sec. |6| this is just given by the difference 
in potential energy between the ridge and the minimum — see (6.13). Since V contains 
exponential factors, evidently this difference is maximal in the limit (prain 0ck-cnd- Thus 
we find 



(/"min/A^Pl Xo 



;9imin/A/pi 



PI 

Xo 



1 + 



Meanwhile Knin is to a good approximation given by K 



I ' mm I ~ 



ck- 



(/"min/A/pi 



Now from ( 7.4[ ) and (7.9) we see that the limiting case 



'ck— end 



(7.14) 



is achieved for 



Xo = p^/'^\/p'M-pi/2. Thus the relevant ratio for condition (6.31) can be written as 



AF 



Xo 



2Mi 



PI 



(7.15) 



which has to be much less than unity. But this is clearly the case since (7.3) implies 



Xo ^ Afpi- It follows that the energy density in x oscillations is parametrically small 
compared to the ghost condensate energy density, thereby ensuring that the NEC is violated 
for sufficiently long to cause a non-singular bounce. 



8 Conclusion 

The many aspects of ekpyrotic cosmology discussed in the paper all revolve around the issue 
of initial conditions. Unlike the hot and chaotic beginning of inflation, ekpyrosis proposes 
a cold and nearly vacuous initial state. Put more concretely, for GUT-scale reheating tem- 
perature, the proto-inflationary patch is 10~^^ cm big, whereas the initial Hubble radius in 
ekpyrosis is of order 1 m. In the absence of a concrete theory of initial conditions, however, 
both represent equally plausible initial states. 
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We have shown that ekpyrotic theory is equally successful as inflation in solving the 
flatness and homogeneity problems of standard big bang cosmology. Indeed, the model allows 
for initial curvature and anisotropic components comparable to that of the ekpyrotic field. 
Starting from these generic initial conditions, the universe emerges in the hot, expanding 
phase with a high degree of spatial flatness, homogeneity and isotropy. Note that this is 
very similar to how inflation addresses these problems — one envisions a universe that 
is essentially smooth and flat over the proto-inflationary patch, thereby allowing cosmic 
acceleration to take over. 

The two-field ekpyrotic model displays a tachyonic instabihty along x-, the direction 

orthogonal to the field trajectory. Since fiuctuations in x by definition coincide with entropy 
perturbations, this instability is essential in generating a scale-invariant spectrum. We cannot 
get rid of it without spoiling the perturbation spectrum. We have shown that the desired 
initial conditions at the onset of the ekpyrotic phase can be naturally achieved with a pre- 
ekpyrotic stabilizing phase. By adding a small mass for x ^"^^ couplings to light fermions, 
X gets stabilized well-before the onset of the ekpyrotic phase for a broad range of initial 
conditions. The mass term is only important at early times — it becomes negligibly small 
during the ekpyrotic phase and therefore does not interfere with the generation of density 
perturbations. 

The analysis of the tachyonic instability in terms of new field variables and x l^ad us to 
propose a simplified and more general version of the scenario. Instead of starting from two 
scalar fields each with their own exponential potential, we can think of one field rolling 
down an ekpyrotic direction whose evolution determines the tachyonic mass of a second field 
X- This occurs in such a way that fiuctuations in the latter acquire a scale- invariant spectrum. 
The (0, x) language allows for far greater freedom in specifying the potential. Whereas the 
earlier (0i, ^2) framework required two steep exponential functions in the potential, here we 
only need one steep potential along (p and one coefficient in a Taylor-expansion — the mass 
term for x — to be such that the curvature of the potential is approximately the same in 
either direction: V^-^^ ~ V^^^. There is, therefore, enormous freedom in specifying the global 
shape of the potential. 

The simplified picture in terms of and x has important consequences for the predictions 
of the model, in particular for the spectral index and the level of non-Gaussianity. The former 
now depends on an additional parameter 8 that characterizes the deviation from V^-^-^ — V^tj>(t>- 
With vanishing 5, a pure exponential for leads to a slightly blue spectrum, in disagreement 
with recent data. Allowing for non-zero 5, however, can make the spectral tilt red even in 
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this case. Turning to non-Gaussianity, we have argued that self-interactions in x have a 
natural cut-off scale of order A = ^/eM■p\. Assuming coefficients of order unity, the resulting 
level of non-Gaussianity is therefore high, corresponding to /nl ~ 1/e. The consequences 
for current and near-future microwave background experiments will be discussed elsewhere. 

These achievements, we believe, together put the New Ekpyrotic scenario on a stronger 
footing as a genuine alternative theory of early universe cosmology. 
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